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Abstract 

An experimental comparison of several operational phase 
concepts is presented. In particular, it is shown that statis- 
tically motivated evaluation of experimental data may lead 
to a significant improvement in phase fitting upon the con- 
ventional Noh, Fougeres and Mandel procedure. The anal- 
ysis is extended to the asymptotic limit of large intensities, 
where a strong evidence in favor of multi-dimensional esti- 
mation procedures has been found. 

1 Introduction 

"The essence of quantum theory is its ability to predict 
probabilities for the outcomes of tests, following specified 
preparations" R] . From a pragmatic point of view the quan- 
tum state represents just our information on the system cor- 
responding to a particular preparation by a classical appa- 
ratus. According to quantum theory this seems to be the 
most complete information. However, the accessibility of 
this information is questionable. Not knowing the prepara- 
tion procedure, one does not know the quantum state of the 
system. There is no way to measure it for a single realiza- 
tion of a quantum system. The situation gets better if an 
ensemble of systems prepared in the same quantum state 
are available. Then it is possible to measure complemen- 
tary observables in different experiments, and the quantum 
state of the system can be inferred. Since real ensembles 
are always finite, only the particular numbers of occurrence 
of different results can be measured instead of probabilities. 
This is a paradigm for an arbitrary measurement. How- 
ever, this scheme could seem purposless unless theoretical 
predictions are compared with experiments. In quantum do- 
main this is not in general easy at all, and in practice many 
sophisticated theories cannot be demonstrated on their ex- 
perimental counterparts. 

The estimation of phase differences in interferometry ap- 
pears to be a nice example of the above mentioned scheme, 
where the predictions of quantum theory can be followed 
by an experimental realization. Optical measurements in 
the domain of classical wave optics are well established and 
belong to the most precise measurement schemes currently 
available. Significantly, such schemes may be analysed in 
the framework of quantum phase. 

Quantization based on the correspondence principle leads 
to the formulation of operational quantum phase concepts 



[y, y. Further generalization may be given in the frame- 
work of quantum estimation theory; the prediction may be 
improved using the maximum likelihood estimation. This 
improvement was tested experimentally in matter wave op- 
tics with neutrons, and a statistically significant improve- 
ment was observed [y. This is a remarkable result, since 
the phase estimation is rather uncertain for neutrons due 
to technical limitations of neutron interferometry, where for 
example the visibility of interference fringes is far below the 
ultimate value of 100%. In the present paper the same theo- 
retical background of optimal phase estimation will be used 
for testing of phase resolution with photons. 

Optical measurements offer many advantages. Current 
optical technology enables us to achieve visibility of inter- 
ference fringes close to unity, and to very precisely set the 
intensities of light pulses at levels deep below one photon 
on average. As the main objective, different strategies for 
accurate phase estimation will be specified and their conse- 
quences for achieved precision will be derived. 

The paper is organized as follows. The mathematical 
tools are reviewed in the next section, where the operational 
phase concepts are naturally embedded in the quantum es- 
timation theory. The experimental setup is described in the 
third section. A comparison of several phase estimation pro- 
cedures based on the experimentally measured data is given 
in the fourth section. Finally, the fifth section deals with 
the phase estimation in the asymptotic regime. 

2 Phase estimation 

The operational phase concepts can naturally be embedded 
in the general scheme of quantum estimation theory p, o 
as was done in Ref. |§, 0, i|- Let us consider the 8-port 
homodyne detection scheme [y, y with four output chan- 
nels numbered by indices 3,4,5,6, where the actual values 
of intensities are registered in each run. Assume that these 
values fluctuate in accordance with some statistics. The 
mean intensities are modulated by a phase parameter 6 
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where A^ is the total intensity and V is the visibility of in- 
terference fringes. The true phase shift inside the interfer- 
ometer 0, which is a nonfiuctuating parameter controlled by 



the experimentalist, should carefully be distinguished from 
the estimated phase shift, which is a random quantity. Here- 
after, the latter is denoted by 9 . This device, operating with 
Gaussian signals, represents nothing but a classical wave 
picture of the original 8-port homodyne detection scheme. 
Equivalently, it also corresponds to a Mach-Zehnder inter- 
ferometer, when the measurement is performed with zero 
and 7r/2 auxiliary phase shifters. In this case, data is not 
obtained simultaneously, but it is collected during repeated 
experiments. Provided that a particular combination of out- 
puts {ns, 714, ?i5, ne} has been registered, the phase shift can 
be inferred. The point estimators of phase corresponding to 
the maximum-likelihood (ML) estimation will be used here 



10^ |11| . In accordance with the ML approach |1^, the 
sought-after phase shift is given by the value, which max- 
imizes the likelihood function. Provided that the noise is 
Gaussian, the likelihood function corresponding to the de- 
tection of given data reads 
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Here the variance a^ represents the phase insensitive noise 
of each channel. A notation analogous to the definition of 
phase by Noh, Fougeres and Mandel p] may be introduced 
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The likelihood function (0) may be maximized by the choice 
of parameters for phase shift and visibility, respectively M 
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Hence the operational phase concept of Noh, Fougeres and 
Mandel coincides with the ML estimation for waves repre- 
sented by continuous Gaussian signal with phase indepen- 
dent and symmetrical noises. These rather strict assump- 
tions are incompatible with the nature of signals encoun- 
tered in experiments; such restrictions would be, however, 
natural in the classical theory. 

The optimum prediction is different for Poissonian statis- 
tics. ML estimation based on the Poissonian likelihood func- 
tion 
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gives optimum values for phase shift and visibility 
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provided the estimated visibility (S) is smaller than unity. 
In the opposite case it is necessary to maximize the likeli- 
hood function (M) on the boundary {V = 1) of the physically 
allowed region of the parameter space numerically. Rela- 
tions (&0) provide a correction of the Gaussian theory with 
respect to the discrete signals. Besides the phase shift, vis- 
ibility of interference fringes and the total input energy can 
be evaluated simultaneously. 

The apparent difference between relations (^-§) and (Q- 
H) represents the theoretical background of the presented 
treatment. Obviously, both predictions will coincide pro- 
vided that there is almost no information available in the 
low field limit A^ ^- 0. Similarly in the strong field limit 
TV — > 00, the phase of the light is well defined and both in- 
ferred values of the phase approach the same value. Possible 
deviations may appear in the intermediate regime iV « 1. 
The test of the difference between (||) and (||) is proposed 
as controlled phase measurement. The phase difference was 
adjusted to a certain value and estimated independently us- 
ing both methods (^ and (||) in repeated experiments. 

To compare two or more phase estimators, some measure 
of the estimation error is needed. Dispersion defined as 
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can well do the job. Here the average is taken over poste- 
rior phase distribution of the corresponding phase estima- 
tor. The dispersion ( pi]| ) is a compact space analogy of the 
averaged quadratic cost function (variance), frequently used 
in the estimation theory |^]. 

The evaluation of the average quadratic cost (|o|) is not 
the only way to compare efficiencies of different estimation 
procedures. Another possibility is to use the rectangular 
cost function 
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This choice of the cost function corresponds to the follow- 
ing evaluation the experimental data. Each sample of data 
consisting of numbers 713,714,715,716 of counted photons is 
processed using NFM formula (0) issuing phase prediction 
0NFM- The relative frequency fg{A6), which is proportional 
to the average cost of the Gaussian estimator {C{9 — 6)), 
characterizes how many times the estimated phase Onpm 
falls within the chosen phase window A6 (confidence inter- 
val) spanning around the true phase shift. The same proce- 
dure is repeated for phase predictions based on the Poisso- 
nian phase estimator (|^) yielding the relative frequency of 
"hits" fp{Ae). The quantity 
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represents the difference in efficiency of the ML and NFM 
phase estimations for the given phase window A6 and given 
input energy N . If this quantity is significantly positive, 
the ML estimation is better than its NFM counterpart. On 
the other hand, if AE is close to zero, both data evaluation 
procedures are statistically equivalent and no discrimination 
is possible. 




Figure 1: Scheme of the laboratory setup. FC - input fiber 
coupler, PCx - polarization controllers, PMx - phase mod- 
ulators, ATT - attenuator, SVRC - output variable ratio 
coupler, Dx - detectors. 



3 Experimental setup 

The laboratory setup (see Fig. |l|) is based on a single-mode- 
fiber Mach-Zehnder interferometer carefully balanced and 
adjusted for maximum visibility. A semiconductor laser 
source (SHARP LT015) produces 4-ns-long pulses with a 
repetition rate of 130 kHz. Initial pulse intensity is about 
10^ photons per pulse. This is decreased by 11 dB due 
to losses in the fibers and other components of the setup, 
and precisely adjusted by artificial attenuation in the pro- 
grammable attenuator (JDS Fitel HA9) to reach the re- 
quired level at the detectors. Input coupler FC (SIFAM) 
divides the pulses between the arms of the interferometer 
(each 4 m long). Both arms contain planar phase modu- 
lators PM1,2 (UTP). Only PMl has been used for phase 
settings, the other is included just for symmetry reasons. 
Both modulators also work as linear polarizers (extinction 
ratio 1 : 10^) to improve the degree of polarization. In- 
put polarization to the modulators is set by polarization 
controllers PCI and PC2. Attenuator ATT in the upper 
arm of the interferometer helps balance the losses in both 
arms of the interferometer to reach maximum visibility. The 
length of the arms is balanced by a variable air gap (AG). 
Polarization controller PCS is used to match the polariza- 
tion in the arms at the output variable ratio coupler SVRC 
(SIFAM). The result of interference is detected using sil- 
icon photon-counting detectors (EG&G SPCM-AQ) with 
less than 70 dark counts per second and quantum efficiency 
of 55%. The signals from the detectors are processed using 
detection electronics based on time-to- amplitude converters 
and single-channel analyzers (EG&G Ortec) and recorded 
by a computer, which also controls the driving voltage of 
the phase modulator, programmable attenuator setting and 
laser operation as well. In this setup we have reached inter- 
ference visibility of up to 99.8%. 

The whole interferometer is placed in a polystyrene box 
to minimize thermal drift of the fringes. After initial warm- 
up, the phase stability of the device is better than 7r/3000 
per second. During the measurement, active stabilization 
of the interference pattern is performed each 5-10 s. 
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Figure 2: The experimentally observed difference between 
dispersions of the NFM and ML estimators as a function of 
the input mean number of photons N for fixed true phase 
9 — 7r/3. Error bars corresponding to 68% confidence inter- 
vals are shown. 

4 Measured data evaluation 

Unfortunately, commercially available photodetectors for 
measurement of weak quantum signals fail to discriminate 
the number of detected photons. Only the presence or ab- 
sence of the signal can usually be detected. The impossibil- 
ity to count photons is circumvented as follows. According 
to the well-known polynomial theorem, the sum of two or 
more Poissonian signals is a Poissonian signal again, the 
mean simply being replaced by the sum of the means of its 
constituents. It is therefore possible to carry out measure- 
ments with very weak signals of intensity, say, 0.01 — 0.001 
photons per pulse so that probability of two photons be- 
ing in the same pulse (double-detection) is very small, and 
then collect an appropriate number of individual yes-no de- 
tections to obtain desired "input" intensity N. For exam- 
ple, an experimental run with input pulse mean intensity 
TV = 10 can be simulated by a sequence of 10000 measure- 
ment with mean input intensity Np = 0.001 photons per 
pulse. The probability of double-detection in a single run is 
p < 10~^ for a Poissonian light source. Hence the probabil- 
ity of single double-detection during the whole sequence of 
measurements is less than 1% and the probability of triple- 
detection or several double-detections is entirely negligible. 
This procedure enables us to effectively simulate the results 
of experiments with intense pulses N ^ 1 and ideal pho- 
todetectors. Whenever in the text an experimental sample 
is mentioned, it should be clear that we actually refer to a 
sum of many experimental samples measured with intensi- 
ties well below a single photon per pulse. 

The difference of dispersions (Eot) of the Gaussian and 
Poissonian phase estimators found in our experiment is 
shown in Fig. ^ for a fixed true phase 6 — tt/3. The number 
of experimental samples used for calculation of the disper- 
sions varies from 1000 samples for input intensity A^ = 60 
to more than 100,000 samples with N = 0.1. The error 
bars arising from a limited number of samples are the result 
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Figure 3: Experimentally obtained AE (squares) compared 
to theoretical values (circles). Error bars corresponding to 
7500 measured samples are shown. 



of numerical simulation. Fig. agrees well with qualitative 
reasoning of the previous section. The most distinct differ- 
ence between the dispersions of the ML and NFM estima- 
tors is seen for the input mean number of photons N w 7.5. 
Thus, it can be said that as long as interference and phase 
measurements are concerned, discrete signals with Poisso- 
nian statistics are distinguishable from the classical wave 
for only a relatively small range of input energies. 

The difference in efficiency of the ML and NFM phase 
estimation (n2[) calculated from experimental data is shown 
in Fig. H The difference AE was calculated using 7500 ex- 
perimental samples measured in experiment with iV = 10 
photons and visibility of 99.6%. The chosen input energy 
roughly corresponds to the maximum seen in Fig. y. Since 
the experimental data are limited to a finite number of sam- 
ples due to experimental conditions and available time, the 
estimated AE would be slightly different in repeated exper- 
iments. Statistical significance of the experimental results is 
demonstrated using computer simulation again. Standard 
deviation corresponding to 7500 measured samples is shown 
in Fig. ^ as error bars for each phase window. 

A significant difference between the effectiveness of clas- 
sical and optimal treatments is apparent in the Fig. pi The 
optimal treatment provides an improvement in estimation 
procedure, and the difference is beyond the statistical error 
by more than 10 standard deviations in the optimum case. 
High stability and visibility of interference fringes in the 
optical interferometer along with a high repetition rate of 
the pulsed laser make the improvement of the NFM phase 
prediction more evident than in a similar comparison per- 
formed with the neutron interferometer setup B. Notice 
the dependence of the precision gain AE on the width of 
the chosen phase window. Obviously, no better performance 
of the ML method can be expected for large values of the 
phase window A0; any sensible statistical method would 
yield quite reasonable results. Likewise, no real improve- 
ment over the Gaussian estimate can be expected when A9 
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Figure 4: Observed relative difference between dispersions 
of the NFM and ML estimators. All the parameters are the 
same as in Fig. g. 

is close to zero, because too few data would then fall within 
the window. The largest difference is about 6% in the win- 
dow of the width of about 0.5 rad. 

5 Asymptotic behavior 

As intensity of the input light increases, both estimations of 
phase shift yield more sharp and precise results. The error 
of the best-known proposed phase measurements scales as 
N~^ for large input intensities ||l^, |lj]. Our photodetec- 
tion scheme cannot compete with such measurements. On 
the other hand these methods necessitate the use of exotic 
states such as two-mode Fock states, etc., which are still 
impossible to prepare in contemporary laboratories. Since 
the experimental equipment is always limited, the only way 
to improve precision of phase measurements lies in careful 
evaluation of the measured phase sensitive data and distil- 
lation of all available phase information. Therefore, it is 
worthwhile to compare the performance of the NFM and 
ML estimations in the limit of high intensities. 

To get some qualitative feeling of how both estimators 
approach the above mentioned limit we redrew Fig. with 
differently scaled vertical axes, see Fig. 0. It can be seen that 
the relative difference of both dispersions monotonically in- 
creases with N and finally approaches some constant value 
different from zero. This means that during the transition 
from TV « 10 to higher values of N, both estimators first 
scale with slightly different powers of N, and for high in- 
tensities both powers reach the same value and the ratio 
of the Gaussian and Poissonian dispersions approach some 
constant. 

It is easy to calculate the behavior of the dispersion of 
the Gaussian phase estimator for N -^ oo. The dispersion 
becomes 
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Readers interested in details are referred to Appendix ^. 
As can be expected, the error of NFM phase measurement 



is proportional to N^^'^. This precision represents the so- 
called standard quantum limit. Unfortunately, it is impos- 
sible to derive a simple expression similar to Eq. (O) for 
the Poissonian estimator. However, such a formula is easily 
obtained provided the physical constraint F < 1 is released. 
Thus throwing the estimated value of visibility away and 
interpreting Eq. (g|) as an estimator of the unknown phase 
shift valid for each sample {713,714,715,716}, the asymptotic 
dispersion of such an unconstrained estimation reads 
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It is obvious that omitting useful information gained from 
the data makes the unconstrained estimation somewhat less 
efficient than the original constrained one. 

We can see from Eqs. (O) and (O) that in the limit 
of low visibility both the NFM and ML phase predictions 
are equivalent. This result is in agreement with properties 
of the well-known discrete Fourier transform (DFT) phase 
estimator ||l^ in the same regime[|. On the other hand, 
ML estimation always gives better results than NFM the- 
ory, provided the visibility is high. For some values of the 
true phase shift, reduction in the dispersion down to 50% is 
possible. 

The ultimate limit to the resolution of the particular es- 
timator is set by the well known Cramer-Rao inequality. 
Provided the visibility and input intensity are under con- 
trol in the experiment, the phase shift 9 remains the only 
parameter to be estimated. For such a single-parameter 
problem, the Cramer-Rao lower bound (CRLB) on the dis- 
persionn of any phase estimator is given as follows 
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Here the symbol E denotes averaging over observed data. 
Upon substitution of the joint Poissonian distribution of 
the sample {713,714,715,716} to Eq. (pq) and making similar 
approximations to those used in derivation of Eq. dl3), we 
end up with 
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Two interesting observations follow from (|l6|). First, notice 
that for low visibility, the dispersion of the NFM phase pre- 
diction (|3|) attains the CRLB. This means the NFM estima- 
tion is best possible in this hmit. Second, for perfect exper- 
imental setup {V=l), the CRLB is simply cr'^m^B = 1/2N. 
Thus the resolution of the unconstrained ML estimation 
can still be improved a bit. In deriving Eq. (Ilq) we sup- 
posed the value of visibility and input intensity of the laser 
beam are known. Such knowledge represents some addi- 
tional information about the experimental setup. Let us 

^We note in passing that DFT can in fact be regarded as a general- 
ization of NFM phase concept to a greater number of auxiliary phase 
shifts [y. 

^Actually, the CRLB holds for variance rather than for dispersion. 
Since all the relevant phase uncertainties are small in the limit of high 
input energy, both quantities coincide in this case. 
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Figure 5: Bias of the single-parameter phase estimator for 
expected visibility V—1 as a function of the true phase shift 
9. Actual visibility of interference fringes varies from V=l 
(straight line) to y=0.1 (most bent line) in steps AF=0.1. 



Estimator 


a' 


Ce 


Eja''d9 


NFM 


1/N 




2tt/N 


unconstr. ML 


(1 + cos2 29)/2N 




|7r/iV 


constr. ML 


« (l + 0.5cos2 26l)/2iV 




«|7r/iV 


CRLB 


l/2iV 




vr/iV 



Table 1 : Asymptotic dispersion and overall quadratic cost of 
various phase estimators. For comparison, CRLB is shown. 
Note that phase prediction of the ML estimation with phys- 
ical constraint on the inferred value of visibility is superior 
to the prediction without the constraint. 



assume that the input intensity and visibility are really un- 
der control and, in addition, the visibility is equal to unity. 
In this case the Poissonian likelihood function (0) depends 
only on the value of the phase shift 9. Now the single- 
parameter ML estimation of the phase shift 9 consists of 
maximizing the likelihood function C{9, V=l) with respect 
to the single parameter 9. This procedure is just what we 
have done in the case of many-parameter ML estimation 
(§r0), when the experimental sample yielded an unphysical 
value of visibility V > \. The only difference is that in 
the case of single-parameter ML estimation we maximize 
the likelihood function on the boundary for any detected 
sample {713,714,715,716}. We may ask whether the single- 
parameter phase estimator achieves the best phase reso- 
lution aQj^j^g — l/2N. An explicit calculation (for details 
see Appendix H) shows this is really the case. Although 
it may seem that the single-parameter ML estimation thus 
gives best results, some caution is necessary when visibility 
(or another parameter) is not known precisely or fluctuates. 
For example, estimation on the boundary V = 1 leads to 
a strongly biased phase prediction provided the actual vis- 
ibility differs from unity, as is demonstrated in Fig. g[ In 
this particular case, the bias caused by dismissing the pos- 
sibility V^ < 1 is independent of the input light intensity N. 
For larger intensities it dominates the uncertainty of esti- 
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Figure 6: Asymptotic dispersion of the NFM estima- 
tor; theory (soUd hne) and experimentally obtained values 
(squares). Asymptotic dispersion of the unconstrained ML 
estimator; theory (dashed line) and experimentally obtained 
values (triangles). Experimentally obtained dispersion of 
the constrained ML estimation (circles) . The corresponding 
input mean number of photons and the estimated visibility 
are N = 160 and V — 99.2%, respectively. 



mated phase and the single-parameter ML estimation may 
be outperformed by the Gaussian (NFM) one. Therefore 
one should always estimate all parameters, which are not 
under experimentalist's control together with the parame- 
ter of interest regardless of smaller theoretical effectiveness 
of such a complex estimation procedure. 

Now let us return to the many-parameter constrained ML 
estimation. In this case a particular detection {713, n^, n^, 
rie} is processed either by Eq. (g) (when applied to all sam- 
ples, (T^ is given by Eq. (M)) or via maximization of like- 
lihood function (H) on the boundary (when applied to all 
samples, a'^=l/2N). Although we do not switch between 
these methods at random, because the choice depends on 
the particular sample, the mean of both dispersions gives 
us a rough estimate of the performance of the constrained 
ML phase prediction. A more precise value can always be 
obtained with the help of computer simulation. The per- 
formances of various phase estimators are summarized in 
Tab. |]for a perfect experimental setup V—1. 

An experimental comparison of the three phase estima- 
tions in the asymptotic regime is shown in Fig. pi For 
comparison, the theoretical values of dispersions given by 
Eqs. (0) and dlj) are also shown. The dispersions were 
determined using 10, 000 measured samples with N = 160 
for each value of the true phase shift 9 — /ctt/IG rad, 
k — 0,1, ..,8. More than 10^ weak laser pulses were sent 
through the interferometer to obtain the figure. Several im- 
portant conclusions can be drawn from Fig. |6[ (i) We can 
see that the uncertainty of the constrained ML estimation 
is definitely below the uncertainty of the unconstrained es- 
timation in agreement with our arguments presented in this 



section. It means that insisting on the physical constraintsn 
of allowed results of estimation or reconstruction procedure 
is important not only for interpretation reasons, but it also 
makes the estimation more efficient, (ii) The observed val- 
ues of dispersion exhibit a systematic error. The additional 
noise above the theoretical uncertainty is caused by inherent 
phase fluctuations in the experimental setup and their mag- 
nitude can be estimated from Fig. g as 0.020 ± 0.003 rad. 
This value is in an excellent agreement with the value 0.019 
rad obtained by an independent method. Hence our sta- 
tistically motivated evaluation of experimental data can be 
used for inferring the amount of fluctuations, and therefore 
it provides an independent and nontrivial way for calibrat- 
ing an interferometer. Moreover, a slightly different sen- 
sitivity of different phase estimators to various parameters 
of the setup makes it possible, at least in principle, to dis- 
tinguish between different sources of noise. This is another 
interesting feature of the method we propose, (iii) In Sec- 
tion |4| we could see that the most distinct difference between 
semiclassical and fully quantum phase concepts occurs in 
the regime, where the intrinsic phase uncertainty of light 
is much larger than phase fluctuations caused by any rea- 
sonable imperfections of the experimental setup. Therefore, 
though clearly visible in Fig. O, the "external" phase fluc- 
tuations may be completely neglected in Fig. |[ However, 
with increasing intensity the (unavoidable) fluctuations be- 
come comparable with the intrinsic phase uncertainty, and 
for even larger N the accuracy of any phase measurement is 
governed by the external influences rather than by the the- 
oretical limit of the corresponding phase estimation. The 
statistics of light are then no longer reflected in its phase 
properties, and different quantum phase concepts become 
indistinguishable. Not only does this provide another evi- 
dence for the fact that the NFM phase concept differs from 
its ML counterpart only for a narrow range of energies, as 
we already stated in section |4[ but it also shows how the op- 
erationally defined quantum phase approaches its classical 
limit. 



6 Conclusion 

Theoretical and experimental justification of operational 
quantum phase concepts is addressed in this paper. Statis- 
tically motivated evaluation of the interferometric setup has 
been presented. The choice of the optimum phase estimator 
strongly depends on the experimentalist's knowledge about 
the interferometric setup and on the nature of the signal be- 
ing detected. Two important cases - NFM and ML estima- 
tors, resulting from the classical and quantum description 
of the experiment, respectively, have been compared. Dif- 
ferences between both treatments have been measured ex- 
perimentally and have shown to be statistically significant 
in the limited range of input energies. In particular, no dif- 
ference between the NFM and ML phase predictions have 
been observed in the regime of very small number of parti- 
cles, which is usually considered as the domain of quantum 
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physics. This detailed analysis makes it possible quantify 
the amount of the noise associated with the phase. The 
lack of knowledge about the parameters of interferometric 
setup has also been considered. In the asymptotic limit of 
large input energy the intrinsically biased ML estimation 
procedures yield sensible results only provided that all the 
uncertain parameters of the setup are estimated together 
with the unknown phase shift. This can be interpreted in 
the framework of more complex estimation procedures - the 
so called quantum state reconstructions. 
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Asymptotic dispersion of 
Gaussian (NFM) estimator 



the 



To calculate the dispersion of the phase estimator, we need 
to evaluate the expectation of the sine and cosine functions 
of the inferred phase, e.g. 



y^ cos [9{n3, n4, n^, ne)] J]^ P{ni). 



(17) 



j=3 



Here the inferred phase shift 9 conditioned by the detec- 
tion (ri3, ^4, 725, rig) is given by Eqs. (||) and (||) for the 
NFM and ML estimations, respectively. Now, suppose the 
interferometer is fed by a strong pulse with A^ ^ 1. Pro- 
vided the true phase shift 9 ^ kir/A, fc e AA, we also have 
fii ^ 1, i = 3,4,5,6, and the Poissonian photocount dis- 
tribution P{ni) can be approximated by Gaussian with the 
same variance, near its peak: 



Pin,) = 



'(" 



,r/2n 



V27m7 
Hi > 1, Tii — hi <^ fii. 



(18) 
(19) 



Notice that although the distribution now becomes sym- 
metric (i.e., the estimation is unbiased for large N), the 
noise remains phase sensitive even for high input energy in 
contrast to the assumption hidden in the NFM theory (S). 
Using the Gaussian phase formula (0) in (O), we obtain 
the expectation value of the cosine phase function in the 
following form 



(COS! 



"-34 



V^li + ' 



--P{n34)P{n5e)dn34dn5G, (20) 



'56 



where 7134 =713 — 714, n^Q ~ n^ — uq and we used the fact 
that the numbers of counted photons appear in the Gaus- 
sian exponential phase estimate (|3|) only in terms of their 
differences. The probability distribution of the differences 



can be calculated from the photocount distributions follow- 
ing the simple rule 



P{n,,) 

where ij 
we have 



P{ni)P{nj)5{ni — Uj — nij)dnidnj, (21) 



34, 56. For Gaussian probability distributions 



^("34) = 



\/27r(n7+n4) 



2(fl3+fl4) 



(22) 



for example. Since the signal to noise ratio is large in the 
limit of high intensity, it is legitimate to split the counted 
numbers of photons (or their differences) into their mean 
values and small fluctuating parts 



A., 



i = 34,56. 



(23) 



Now we expand the estimated sine and cosine phase func- 
tions keeping only the fluctuation-independent term and 
the second-order terms in the fluctuations. For cos 9 we get 
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2' 



3 ?^34nig 2 



{nl 



A 



\5/2 "56- (^4) 

'34 T "56^ 

The expansion of sin 9 is obtained exchanging 34 <-> 56. 
Substituting Eqs. (H) and (H) to Eq. ^) and using the 
following relations 



(A34) =»^3+n4, {Al^)=n5+ne, 


(25) 


nl^^N'^V'^, 713 -f 774 = 775 + Tig == A^, 


(26) 



"■34 

the first set being implied by Eq. (E2), we arrive at an ap- 
proximate mean value of the cosine function of the estimated 
phase 



1 



2NV^ 



(cos^ 



'■9). (27) 



An analogous expression for sin 9 is obtained exchanging 
cos 9 ^-> sin 9 in Eq. (|2j) . Finally, using Eq. (g7|) and a sim- 
ilar expression for sin 9 in the dispersion formula din) and 
neglecting terms of the order l/N"^, we arrive at the asymp- 
totic dispersion of the Gaussian estimating procedure (^3). 
Since a finite change in the effectiveness of the estimation 
caused by an infinitely small change of the estimated param- 
eter is unphysical, the derived expression for the asymptotic 
dispersion also holds for the isolated values of the true phase 
for which our procedure fails. 

In the case of ML estimation we can proceed in a com- 
pletely analogous way. Starting from the expansion of the 
Poissonian phase estimator (0) in fluctuations of the num- 
ber of counted photons, we obtain, by a straightforward but 
a rather lengthy calculation, desired expectation values of 
the cosine and sine phase functions 



(cos 9)kcos9 + Y^ C, {9, V, N)f 

1=3 
6 

(sin 9) K sin 9 + '^S^{9,V,N)f 



(28) 
(29) 



i=3 



Here Q's and Si's are the coefficients of the terms of the 
Taylor series, quadratic in corresponding fluctuations A^. 
When an explicit form of the coefficients is substituted into 
( pq ) , it is then easy to obtain the asymptotic dispersion (Q) 
of the unconstrained Poissonian ML estimator (ph. 

B Asymptotic dispersion of the si- 
ngle—parameter Poissonian esti- 
mator 

As above, we will suppose that detected numbers of photons 
can be decomposed into their mean values and fluctuating 
parts, small with respect to the means, as follows, ni=ni + 
Arii. Inspection of the Poissonian likelihood function (0) 
shows that the point is a local maximum of C if and only 
if the condition 

f-\n[C{e)]^0 (30) 

holds. Assuming a perfect experimental setup, V—1, the 
derivative of the log-likelihood function becomes 



^'logiO) = -«3 
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\- cos 9 
cos 9 
1 — sin ( 
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sin 6* 
1 — cos 9 
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cos6' 
l + sm9 



(31) 



where C'l^g = d\n{C)/d9. Now we make use of the fact that 
the result of unconstrained ML estimation (||) is not so bad 
and in particular it lies close to the true global maximum of 
the likelihood function (0). The purpose is two-fold. First, 
the estimated phase (||) is a good starting point for finding 
a root of the expression (^) with the help of some approx- 
imation method; second, it automatically selects the global 
maximum of C among all possible roots of Eq. ( pO[ ) . 
In order to improve our initial guess ^o^ 



0o = arg{e*«} 



(32) 



where exp(i6') is given by Eq. (^), we use the Newton 
method. Since the second derivation of log-likelihood al- 
ways differs from zero, the algorithm converges quickly and 
one step of the method is usually enough to find the global 
maximum with sufficient accuracy. The improved value of 
the estimated phase shift thus reads 



C'i 



C"( 



(33) 



Now we expand the sine and cosine phase functions around 
^0 keeping terms up to the second order in the correction. 
The expectation values of the two functions become 
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(34) 



Finally using Eqs. (||), (|), (|32|) and (|l|) in (||), expanding 
the result to the second order in fluctuations A^, replacing 
Poissonian photocount distributions by Gaussian with the 
same variance, carrying out the average, substituting the 
resulting expectation values of the sine and cosine phase 
functions into Eq. (Eoh, expanding the dispersion in 1/A^ 
and keeping terms at most linear in 1/iV, we arrive at the 
asymptotic dispersion of the one-parameter ML phase esti- 
mation in the form 



a^ = 



1 

2iV 



O 



1 

iV2 



(35) 



The enormous amount of calculation work necessary to ob- 
tain this result was carried out with the help of Maple V 
symbolic mathematical language. 
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